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Abstract 



A class of purely- nonperturbative (PNP) composite operators is defined in 
Quantum Chromodynamics, which is perturbative scheme and scale indepen- 
dent and are useful to describe the internal structure of a strong interacting 
system. The operator product expansion in terms of the new operators cleanly 
separates the perturbative and nonperturbative physics without introducing 
any factorization scale. A number of PNP observables of the nucleon is briefly 
discussed including the PNP parton distributions. In particular, the fraction 
of the nucleon momentum carried by the purely-nonperturbative gluons is 
found to be around 16%. 
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Composite operators composed of elementary fields are used frequently in quantum field 
theory [|l|. They are particularly useful in constructing effective field theories or operator 
product expansions (OPE) when multiple scales appear in a physical problem. It appears 
that there are two classes of intrinsically different composite operators: those that are as- 
sociated with the conserved currents arising from continuous symmetries of a quantized 
lagrangian, such as the vector currents in Quantum Chromodynamics (QCD), and those 
that are not. The latter type of operators, in their bare form, are typically divergent and 
require renormalization. The renormalized operators are known to have scheme and scale 
dependence. One argues that these operators are present only in intermediate stages of a 
calculation and the final physical results are independent of any particular renormalization 
procedure. 

Nonetheless, there are strong motivations to go beyond the viewpoint that composite 
operators are just auxiliary quantities. In an OPE, one often views the composite operators 
describing the physics below the perturbative renormalization scale /i^. When /x^ is at the 
hadron mass scale, these operators act as a probe into the nonperturbative wave function 
of a hadron system. In particular, the nucleon matrix elements of composite operators are 
long regarded as useful observables exhibiting the structural information of the nucleon. For 
instance, the scalar charge of the nucleon (P | ?/''?/' |P), quark distributions, etc. are routinely 
calculated in nucleon models. 

However, scheme and scale dependence of composite operators makes the above interpre- 
tations unconvincing. For instance, it is unknown how much physics in the matrix element 
of a renormalized operator Or{^) is attributed to perturbation theory and how much is truly 
nonperturbative, i.e., a reflection of the bound state properties. The conventional wisdom is 
that the perturbative contributions can be evolved away by going to small /x. In practice, no 
consensus has ever been reached as to precisely which. Strong sensitivity of matrix elements 
in low scale renders any particular choice difficult to justify. Furthermore, a renormalized 
composite operator is arbitrary up to a multiplicative perturbation series 1 + YltLi Cn^s (a*^), 
which can be significantly different from unity when becomes large. A less severe, but 

still significant problem is that the experimentally-extracted matrix elements of composite 
operators, such as parton distributions, cannot be seriously compared with model calcula- 
tions. In a nucleon model, the physics associated with renormalization is often obscure; it 
is difficult to assign a renormalization scale and/or scheme to model matrix elements [^. 

In this Letter a new renormalization definition is offered for a class of logarithmically- 
divergent operators in QCD. The newly-defined operators are perturbative scheme and scale 
independent. Their matrix elements in a hadron state depend only on the strong interac- 
tion scale Aqcd and therefore are totally nonperturbative. To distinguish them from the 
conventionally-renormalized operators, we call them "purely-nonperturbative (PNP) oper- 
ators". An interesting by-product is that the coefficient functions in a reformulated OPE 
with the PNP operators are also perturbative scheme independent. 

The Letter is organized as follows. First the new definition of the renormalized composite 
operators in QCD is presented. Then these operators are shown to be perturbative scheme 
and scale independent. Following that, their use in OPE is discussed and the corresponding 
coefficient functions are shown to be also perturbative scheme independent. Subsequently, 
the new form of OPE and the traditional momentum factorization are compared. Finally, a 
number of comments are made about applications of the PNP operators. In particular, the 
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notion of the PNP parton distributions is introduced. 

Let us begin by recalling the standard definition of composite operators in quantum field 
theory 0. For simplicity, we will not consider explicitly the operator mixing, for which the 
discussions below are unchanged except all equations must be in matrix form. Consider 
a bare composite operator Ob, and insert it into all the renormalized elementary Green's 
functions. Assume primitively-divergent Feynman diagrams have only overall logarithmic 
divergences, which have to be subtracted to define the renormalized operators, 

dnin) = Z-\fi)dB , (1) 

where Z{fj,) is a perturbation series in as{f^). From the viewpoint of factorization, Z{fj,) 
factorizes out the contributions to the bare operator from momentum scales larger than n. 
Oii{fi) is far from unique because it not only depends on the scale fi but on choices of the 
finite part of the renormalization constant Z{fi). 

The central point of this paper is that in QCD one can define a finite operator O^^^^ from 
C_r(/u) by making a further factorization: 

6p-p = Az-\f,)dM , (2) 

where AZ(/i) is defined as 

AZ(,) = („,(,))W^* exp (- + ?1 ?1 . (3) 



where 7(0;) = Z~^ fidZ / d/i = jQa/in + ... is the anomalous dimension of the operator and 
/?(«) = — /?ott/4vr + ... {Pq = 11 — 2nf/3) is the QCD beta-function. This additional, 
perturbatively-calculable factorization takes away all the dependence on the scale fi, intro- 
duced artificially in the perturbative factorization in Eq. (1). AZ{fj,)Z{fj,) is /i independent 
as one can check by direct differentiation. O^^"^ can have no other scale dependence except 
Aqcd and therefore is a purely-nonperturbative operator. Notice that generally O^'^^ cannot 
be defined order-by-order in perturbation theory. In particular, its matrix element in any 
perturbative state takes the form. 




7oln^ +C0 + 



(4) 



where rriq is a soft quark mass and Cq some constant. 

OP'^P is also independent of perturbative factorization schemes. Assume that one has 
adopted a different factorization, 

dn = Z-\^?)dB , (5) 

where Z{^'^) differs from Z{fi'^) by a perturbation series z{a) = 1 + ZiO^. The new 
anomalous dimension becomes 

7(a) = 7(a) + A-A^ . (6) 

Substituting 7(a) in AZ, one has 
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Az-\fi)dR{fx) = /\z-\^)dM = 6p°p . (7) 

Thus apart from an overall multiplicative numerical constant, 0^°^ is uniquely defined. If 
one further imposes the condition that for a conserved current, 0^"^ must coincide with 
the bare operator, the multiplicative constant is reduced to a function of 70, which takes 
the value 1 when 70 = 0, e.g., (A)"^"^/^^"^. This freedom is intrinsically different from the 
ambiguity in perturbative factorization schemes, and is analogous to the arbitrariness in 
the choice of a unit for electric charge and the corresponding electromagnetic field. The 
"natural" convention, as we will see below, is the one we have taken for AZ in Eq. (3). 
One can obviously generalize the above discussion to the quark mass renormalization. The 
result is the so-called renormalization group invariant mass, which has been known in the 
literature for some time [^]. 

The measurement of the matrix element of a PNP operator generally requires high mo- 
mentum processes. Consider a physical observable A{Q), where is a large physical scale, 
such as the mass of a virtual photon. Assume A{Q) has a conventional OPE with a leading 
term depending on M{fj,), the matrix element of Or{^), 

A(g)=C(gV/i^a.(/x))M(/x) + ... , (8) 

where C is a coefficient function calculable as a perturbation series in Q;s(/i). [We ignore the 
possibility that the C perturbation series may be divergent due to infrared and ultraviolet 
renormalons or instanton-anti-instanton pairs 0.] Clearly both M(/i) and C are factor- 
ization scheme and scale dependent. However, because A{Q) is a physical observable, the 
coefficient function C obeys the following renormalization group equation, 

dC 

/i^-7(as(/i))C = . (9) 
The solution of Eq. (9) may be written as, 

fasiQ) ^(^a) da\ 



C[Qyfi\a,{fi)) =C(l,a,(Q))exp - £ 



^c(i,.«,)M«r-e.p(-/;-(i|).gg 

In the second line of the above equation, the Q and fi dependence is completely factorized. 
The /i-dependent part is "naturally" the AZ~^ defined in Eq. (3). Therefore we have, 

A{Q) = C'-^\as{Q))M-^-'^ + ... (11) 

where M^^p is the matrix element of Op'^p in Eq. (4). Now C™^ is a function of as{Q) only, 

C'-(a.(Q)) = Co (a.(Q))^»/2* f 1 + £ c„<(Q)') , (12) 



n=l 
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and is perturbative scheme independent although the couphng constant as{Q) is not. This 
means that C™^ is as physical as the perturbation series for the total e'^e~ — > hadron cross 
section. Through this physical probe C™^, the nonperturbative matrix element MP'^p can 
be extracted from the observable A{Q) in the asymptotic limit, 

MP-^P = lim A(g)/C''^^(a,(g)) . (13) 

Equation (10) represents a nonperturbative factorization: all probe-related (Q^) physics 
is included in C"^(as(Q)), which is generally nonperturbative in as{Q) (70 7^ 0), but is per- 
turbatively calculable; all target-related physics is in MP"p, which is purely-nonperturbative, 
ultraviolet safe, and unique. Emergence of the renormalization group invariants from solv- 
ing evolution equations of the coefficient functions has been noted before 0; however, its 
significance, to the author's knowledge, was not explored. 

At this point it is useful to comment on the relation between this new form of factorization 
and the traditional momentum factorization. If one computes A{Q) in perturbation theory 
in a quark-gluon state, mass singularities may arise in the form of the following integral, 

where niq is a quark mass. To isolate the contributions from the small k± region, one usually 
splits the integral by introducing a factorization scale fi, 

,1^ kj_ Jm^g kj_ 

The first term contains hard momentum only and can be identified as a part of the coef- 
ficient function C. The second term is soft and is identified as a part of the perturbative 
matrix element of a nonperturbative operator 0^(/i). It seems that the asymptotic free- 
dom is best respected in this way. According to the discussions earlier, however, the second 
term still contains perturbatively calculable contributions. Indeed, a truly nonperturbative 
quantity, such as hadron mass, should not be sensitive to the ultraviolet cutoff fi. The 
present /i dependence in the second term can be further factorized away, leaving the purely- 
nonperturbative contribution, as signalled by the mass singularity, intact. However, this 
cannot be done order-by-order in perturbation theory as AZ in Eq. (3) and the perturba- 
tive matrix elements of the PNP operators in Eq. (4) are generally non-analytic in as{n). 
An order-by-order isolation of the mass singularity necessarily generates the perturbative 
/i-dependence in the soft part. From a different perspective, one might view the momentum 
factorization as the process of splitting the physical C™^ into a product of two perturbative 
factors C(/i) and AZ{fi), and then combining the latter with the the matrix element in Eq. 
(4) to produce a perturbative expansion in as(/i). To be sure, this might be a necessary 
intermediate step to calculate C"^^. 

It is interesting to note that for a given conventional factorization scheme and a given 
operator, one can always find a /io such that AZ(/io) = 1. At this fiQ, the above momentum 
factorization into perturbative and nonperturbative contributions seems effectively complete. 
Therefore, one could take /io as a matching point between model calculations and data. Since 
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/io is generally small (^^(yUo) ~ 1), the actual perturbative evolution to such low scales is 
questionable. 

The usefulness of the PNP operators lies in their applications to studying the internal 
structure of hadrons. As we have said, the operators with given quantum numbers and 
dimension are essentially unique and completely nonperturbative; therefore, it is proper to 
interpret them as quasi observables. On the practical side, the matrix elements of various 
operators calculated in nucleon models are most naturally identified with those of the PNP 
operators in QCD. [Of course, the issue of "effective charges" still remains 0.] The observ- 
ables measured in high energy processes, such as parton distributions, can now be totally 
decoupled from the high-energy aspects of the probes and attributed solely to the internal 
structure of targets. Before closing, a few comments are made in this direction: 

1) Parton distributions are normally defined in a scheme (usually MS) and scale. One 
can now define purely-nonperturbative parton distributions of hadrons with their moments 
equal to the matrix elements of the PNP twist-two operators. For many years, Brodsky 
and his collaborators have studied the "intrinsic parton distributions" [§. They assert 
that intrinsic distributions can be calculated in terms of the (light-cone) wave functions 
of the hadrons only and have a weak scale-dependence. The intrinsic distributions may be 
identified with the PNP distributions. Compared with the ordinary parton distributions of a 
given scheme, the PNP distributions have a tamed small-x behavior. Indeed, the AZ factor 
contributes to much of the small-x rising in the ordinary parton distributions because 70 is 
a monotonically increasing function of the spin of the twist-two operators and is singular or 
negative near spin-zero. The results for the PNP distributions, together with other studies, 
will be presented in a separate paper 0. 

2) According to standard perturbative calculations, gluons carry a fraction 16/(16 + 3nj) 
of the nucleon momentum and spin in the asymptotic limit |T^. Although this result is 



scheme independent also, it is completely perturbative as it is independent of hadron states. 
The fraction of the momentum carried by the purely-nonperturbative gluons, 71"^"^, is related 

to the conventional fraction 7r^^(/i) hj [uf = 3), 



nr^' = 0.64 + 



_0.64H-0.12^ + ,f(,) l-0.57«^) 



TT ^ \ TT 



a.r(,)) (16) 



at the next-to-leading logarithmic order. Using 7r^^(1.6 GeV) = 0.42 |]TT| and a^°(1.6 GeV) 



= 0.3, one finds that the gluons carry only 16% of the nucleon momentum. This not only 
confirms that the PNP glue distribution is much reduced, but also lends important support 
for valence quark pictures. 

3) The fraction of the nucleon spin contributed by the spin of the quarks AE(yu) (singlet 
axial charge), defined in {P\ J2i=u,d,si'i'l^'l54'i\P) = 2AS(/i)S'^, is scheme and scale depen- 
dent because of the Adler-Bell-Jackiw anomaly. The purely non-perturbative fraction can 
be defined, 

where 7(a) starts at order a^. This definition of EP'^p has already been recognized before 
HJll. The actual deviation of AEP'^p from, say, AE^(/i2 = 10 GeV^) is small. The 
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fraction of the nucleon spin carried by the purely-nonperturbative gluons in the hght-hke 
gauge (A"*" = 0) is related to the conventional scheme and scale dependent fraction by, 



^QPnp _ ^ cia,(/i) + ...)as{fi){AG{fi) + 4AS/6etoo) (18) 



Using AG(1 GeV) = 1.6±0.9 fTl, we get AGp°p = 0.72±0.4 in the leading order, a reduced 
but still significant number. 

4) In recent years, the so-called strange content of the nucleon has been studied exten- 
sively. While the strange radius and anomalous magnetic moment of the nucleon, defined 
from the strange vector current s'-f^s, are scale and scheme independent, other strange 
observables, e.g., strange scalar, axial, and tensor charges, are not Jl^- Because of this, 
questions have often been raised whether the latter quantities can faithfully describe the 
strange content of the nucleon. One can now use the PNP operators to achieve this. For 
instance, the strange scalar charge defined in {P\ss\P) = s(/i^)(2M) has a PNP version, 

sP'^P = + 0.90^^ + . (19) 

TT 

which can be extracted from the tt — N a term. 
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